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Abstract. We prove a uniqueness theorem for compactly supported initial data 
for the variable speed wave equation arising in models of thermoacoustic tomography, 
given measurements on a part of the boundary. The proof is based on domain of 
dependence arguments and D. Tataru's unique continuation theorem. 

1. Introduction 

Thermoacoustic tomography (TAT) is a recently developed technique in medical 
diagnostics. The idea is to combine the high contrast of electromagnetic radiation 
with the high resolution of ultrasound. In utilizing this combination one hopes to 
overcome the weaknesses of both imaging methods. 

In TAT a brief duration microwave pulse irradiates the object. (The practice of 
photoacoustic tomography uses a laser instead of a MW pulse, and then proceeds 
similarly.) Although the irradiation is assumed to be uniform, more energy may 
be absorbed in some locations than in others. For example, cancerous cells absorb 
several times more energy than healthy cells. For this reason, knowing the energy 
absorption f{x) at each point is of importance as a diagnostic tool. The goal of TAT 
is to determine f(x) by using transducers on the surface of the object to measure the 
pressure variation u(x,t) generated by the absorption of energy. Here u(x, t) solves 
the wave equation 

Pu = u tt - c 2 (x)Au = 

u(x,0) = f{x) 
(1) u t (x,0) = 0. 

Mathematically, the problem is to recover f(x) from knowledge of u(x,t) on the 
boundary (or part of the boundary) of the region inside the transducers. 

The acoustic speed c(x) is assumed to be known and smooth, having been prede- 
termined from ultrasound experiments. In the case of constant c there are several 
methods known for reconstructing /, at least in the case of complete boundary data 
(see |KuKy| for a survey of results). When c is variable, less is known. A recon- 
struction formula was obtained in AgKu using eigenfunctions of the Laplacian in 
the region Q inside the transducers. 

We are interested in the case of variable c and limited (partial) boundary data. 
We think it is unlikely that an exact reconstruction formula for / can be obtained 
for general measurement surfaces; however with knowledge of u(x,t) on a sufficiently 
large subset of the boundary one can conclude the uniqueness of /. 
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2. Domains of Dependence 

Let be a bounded region in R™ with smooth boundary. Let c(x) be a positive 
smooth function on R™ such that < c(x) < M for some number M > 1. Let 
/ G C°°(R n ) have support inside 0. 

The function c(x) induces a Riemannian metric on R™ by setting g^ (x) = c~ 2 (x)5ij . 
This metric differs from the Euclidean metric on regions where c(x) is not identically 
equal to 1. (In many models of TAT it is assumed that c(x) = 1 outside a compact 
region; we did not find this assumption to be necessary for our purposes, so we have 
omitted it.) The length of a curve r : [a, b] — > R" in this metric is 



. <*•(()) 

This length depends only on the image of r in R™ (i.e. is independent of parametriza- 
tion). The geodesic distance from x to y is defined to be 

(3) d(x,y) = ML(r) , 

r 

the infimum being taken over all curves from x to y. 

If S is a submanifold of R™ define the induced geodesic distance ds(x, y) for x, y G £ 
by taking the infimum only over curves that lie in S. Note that d(x,y) < ds(x,y). 
We then have the following uniqueness theorem: 

Theorem 2.1. Let f G C°°(R™) have support inside 0. Let T C dfl and suppose T 
satisfies: 

(P) For any x G there exists p G T such that d(x,p) < dRn\n(p, dfl\T). 
Suppose that T is large enough that T > d(x, L) for any x G 0. Then if u satisfies 

(4) u tt -c 2 (a;)Au = on I"x R + 

(5) u\ t=0 = f, u t \ t =o = . 
// in addition 

(6) w|rx[o,T]/ = 
we have f = (and hence u = 0). 

Note that a solution of (4)- (6) extends to all of R™ x R t by making u even in t. We 
then have u|r x [-t,t] = 0. 

The demonstration of this has two steps. First, we will apply Green's Theorem to 
show that u must vanish on the domain of dependence in [0, T] x R"\f2 determined 
by T x [0, T]. Second, Tataru's Uniqueness Theorem will allow us to conclude that u 
vanishes inside ft also, and hence / = 0. 



Proof. : Consider the following exterior initial boundary value problem: 
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Vu-Av = on (R n \il) x [0,T] 
v\t=o = v t \t=o = 

(7) w|[o,T]xan = 9 

where g is smooth (the solution v(x,t) exists for t < T and is unique). Suppose in 
addition that g = on [0, T] X T. If g is taken to be equal to the restriction of u 
solving (4)-(6) to dfl x [0,T], then v = u on (E n \0) x [0,T]. 

Fix a point p G M n \0. We need the following fact about the function d R n\n(x,p): 

Lemma 2.2. d,Rn\Q(x,p) is a Lipschitz function of x onW l \^l, and Vd Rn \f l (x,p)\ < 
c~ 1 (x) a.e. 

This lemma is proved in [FinPaR ak] in the case c = 1, and we extend their proof. 
Proof. It is immediate from the definition in (|3|) that 

(8) \d Rn \ n (x,p) -d Rn \ n (y,p)\ < d Rn \ n (x,y) , 

so it suffices to show that for every q £ R n \f2 there is a small ball B(q, p) (defined with 
respect to the Euclidean metric) such that d R n\Q(x, y) < C\x — y\ for all x, y £ B(q, p). 

First, if q ^ <9fi, let p be small enough that B(q,p) does not intersect dil. Then if 
r is the Euclidean segment from x to y in B(q, p), we have from ([3|) 

(9) L(r)<\x — y\ sup c _1 (z) . 

z£B(q,p) 

If g G dfl, then for p small enough we can assume without loss of generality there 
exists a smooth function <f> — (f)(xx, . . . , x n -\) such that 

(10) VL n B(q,p) = {x = (xi, . . . ,X n ) | x n > cf>(xi, . . .,!„_!)} . 

As before let r be the Euclidean segment from x to y. If r does not intersect fi, the 
same analysis applies, and d Rn \Q(x,y) < C\x — y\. If r has a portion entering fl for 
the first time at a and leaving it for the last time at b, then 

(11) d Rn \ n (x, y) < C\x -a\ + C\b - y\ + d Rn \ n (a, b) < 2C\x - y\ + d R n\ n (a, b) . 

Let 7 be the projection of the segment from a to b onto dfl parametrized by s, < 
s < 1. Then 

7(s) = (1 - s)(ai, . . . ,a„_i,0) + s(6i, . . . ,6„_i,0) 

(12) + (O,...,O,0((l-s)ai + s6i,...,(l-s)a n _i+s&„_i) . 
Hence 

< (n-l)|6-a| + (n-l) sup |V0||6-a| <C*|6-a| . 

We conclude from this that d^n\n(x,p) is a Lipschitz function. By Rademacher's 
Theorem (see |EvGa] ) . d R ™\n(x,p) is differentiable almost everywhere. 



(13) 



f/7 
ds 



(s) 
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To estimate the gradient almost everywhere, let x 6 R ra \f2 be a point where the 
gradient exists. Fix 6 > and let e be so small that B(x, e) C R™\17, and also 
swp^ z _ x ^ <e c~ 1 (z) < c _1 (x) + S. Suppose \x — y\ < e and let r parametrize the 
(Euclidean) segment from x to y by the interval [a, b]. Then 



\dR»\n(x,p) - d Rn \ n (y,p)\ < \d R n\ n (x,y)\ 

i b m dt 



a c(r(t)) 
< \x — y\ sup c~ 1 (z) 

\z—x\<e 

(14) < \ X -y\( C - 1 (x) + S) . 

Since this holds for every <5 > 0, we must have 'Vd Rn \fi(x,p) < c^ 1 (x) a.e. 



□ 



Continuing with the proof of Theorem 12. 1[ let xo 6 ft. By Property (P) there 
exists y £ dft such that d(x ,ya) < d Rn \ dn (y ,T) < T. Let d Rn \ m (y ,T) = H and 
select a point p 6 R"\£! close enough to yo such that the intersection of the set 

(15) U={(x,t) eF\fix [0,T]\d R n\ n (x,p)+t<H} 

with cMl x [0,T] is contained inside T x [Q,T]. Such a p can be chosen precisely 
because of Property (P). We wish to show that v solving (|7|) vanishes in U. Such 
a set U is called a domain of dependence (sec [KaKuLa ). A slight modification of 
the usual Green's Theorem argument will allow us to make this conclusion. The 
argument must be modified because, as the authors point out in [FinPaRakJ, dU is 
only a Lipschitz surface. In order to overcome this difficulty, we will approximate U 
from the inside by regions with piecewise smooth boundary. 

d R n\n{x,p) is a Lipschitz function, so it can be approximated by a smooth func- 
tion in the sense that, for any e > there exists a smooth function / £ on R ra \f2 n 
{d®n\p(x,p) < T} such that \d&n\ a (x,p) - f e (x)\ < e and |V/ e (x)| < c~ 1 {x) + e (see 
|Az) ) . We can also approximate (1 — 6)d Rn \n(x,p) by a smooth function f ei $ so that 
\(l-S)d M n\ n (x,p)-f etS (x)\ < e and |V/ e ,j(x)| < (1 - S)^ 1 ^) + e for any S > 0. For 
e sufficiently small relative to 8, we have the important estimate 



(16) |V/ e ,ff(s)| <c- x (s) . 

We next define the domains of dependence on which we will apply Green's Theorem. 
For < h < H let R = R eA h = {(x,t) G R"\^ x [0, h] \ f e ,s(x) +t < h}. Then R is 
a domain with piecewise smooth boundary, and for e sufficiently small, we still have 
OR fl d£l C r. In fact, OR consists of three pieces (Figure 1): 

• Si: A slice of {t = 0} 

• S2: A spacelike surface which is part of {(x, t) | f e ,&(x) + 1 = h] 

• S 3 : A slice of V x [0,h] 

Here, spacelike means c(x)|iVr| < |iV t | for all normal vectors (N x ,N t ). The proper 
terminology should be .spacelike or null, but we shall simply write spacelike for brevity. 
Using Green's Theorem (see for example |Tay| ) we have 



= / ut(c~^(a;)u« - Aw) dV dt 
Jr 

= J c ~ 2 ( x )§- t (l v t^J dVdt + J (W x v t ,V x v)dVdt 

- / div x (v t S7 x v) dV dt 
Jr 

(17) = If {c- 2 {x)v 2 t +{V x v,V x v))uj- [ v t p-dStdt. 

1 iSiUS 2 JE2US3 0v x 

If dS* is the measure on dR induced by dV dt, then u> = N t dS and dSt dt = \N X \ dS 
where (N x , N t ) is the outward unit normal to dR. As a result, 

(18) £ (^c- 2 {x)v 2 t +\V x v\ 2 )\N t \-2v t ^-\N x \^ dS = . 
v, and hence Vt, vanishes on S3, so 

(19) / v t —\N x \ dS = . 

Js 3 ov x 

On the other hand, £2 is spacelike, so since we have 



(20) 



c (x)v t T— 
av x 



< c~ 2 {x)v 2 + \V x v\ 2 , 



we obtain a positive-definite form in (yt, V^i;): 



(21) (c- 2 (x)v 2 + |V^| 2 )|7V t | - 2v t ^\N x \ > 

dv x 

with equality only when (y t , V^u) = 0. Therefore (jTSJ) implies 



(22) 



jf ^ C - 2 (xK 2 + |V^| 2 )|7V t |-2« t |^|iV :c ^ dS = 0. 



As a result, («t,V x ii) = on £2. Letting ft, vary from to H, observing that the 
surfaces £2 (ft) sweep out R ei s,H, and using that the initial data vanish, we conclude 
that v = on R e ,d,H- 

Now let (x, t) 6 U. Then dRn\n(x,p) + 1 < H, so for <5 sufficiently small, 
(1 — S)dftn\n(x > p) +t<H. As a result, 

(23) d a . n \n(x,R e ,s,H) < e 

so by the continuity of i>, v = on [/. 

Extending v to be even in time, we also have v = on — i? := {(t, x)|(— f, a;) €E i?}. 
If we shrink H by any sufficiently small amount, then there exists p > such that 
t> = on B(p,p) x [—H,H] while maintaining d(xQ,p) < H. This is the fact we will 
require in the next section. 
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3. Tataru's Theorem and Unique Continuation 

In this section we will make use of the following special case of a uniqueness theorem 
of D. Tataru: 

Theorem 3.1. [Tataru] Let u G D(K, H 1 (W n )) be a solution to 

n 

(24) d 2 u- ^2 aij(x)d Xi d Xj u = 

*>j'<i 

where the coefficients are smooth. Suppose S is a noncharacteristic hypersurface 
containing a point y, and that u vanishes on one side of S near y. Then u must 
vanish in a neighborhood of y. 

We will use Tataru's Theorem to prove the following fact about solutions which 
vanish in a cylinder: 

Lemma 3.2. Suppose u solves Pu = and u vanishes on the cylinder B(z,p) x 
\—D,D]. Here, B(z,p) is defined with respect to the metric g, i.e. B(z,p) = {x G 
R" | d(x, z) < p) . Suppose D is less than the injectivity radius of the metric gij at z. 
Then u vanishes on X = {(x, t) \ d(x, z) + \t\ < D} (see Figure 2). 

Proof. Let \ — Xe G Co°(K) be supported in [— e, e], satisfy x(t) — 1 f° r \t\ < § i have 
X{t) < \ f° r M > | e an d have |x'(£)| < f for all t. Define for < r < D and e small: 

(25) X r , e = {(*, t) | d(x, z)- + - X (t)) + \x(t) < D} . 

Note that for r < p, X r e C B(z, p) x [-D, D], so u vanishes on X r ^ for r < p. 

dX r:<L is a smooth hypersurface except at the points (z, ±D). Let (x, t) G dX r ^, \t\ < 
D. We can compute the normal to dX r ^ f at (x, t) with respect to the metric g + dt 2 
by considering the gradient of the function 

(26) F(x, t) = d(x, z)j + |*|(1 - X (t)) + |x(*) . 

Since x lies inside a system of normal coordinates about z, we have \V x F(x 1 t)\ = 
V x d{x,z)^r — c _1 (x)-^. We can also estimate d t F(x,t) by computing 

(27) d t F(x,t) = (l- X (t))^- X '(t)(\t\-^) . 
If \t\ > e, \d t F\ = 1. If |e < \t\ < e, then 1 - x (t) > \, and 



(28) 



\d t F(x,t)\ < 



i - x(t) - 



3 e 
7 ' 2 



< 1;. 



If |e < |*| < |e, then |*| - |e < \e, and 



(29) 



|d t F(x,t)|<| X '(t)|(|t|-^) <-'7<!- 



e 4 



Obviously if |*| < |e, d t F(x,t) = 0. 
The symbol of the operator P is 



(30) p(^e,r)=r 2 - C 2 (x)|C| 2 
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At the point (x, t), the normal to dX r ^ e is (x, t, V X F, d t F) and satisfies 



(31) p{x,t,V x F,d t F) < l-c 2 (x) [c- 2 {x) — J <0 

so 9X r!e is noncharacteristic at (x,t). 

Let i? = sup,, {r | u|x r , e = Ve > 0}. By continuity u = on Xr i<: Ve > 0. But if 
R < D, the boundary of X^.t is noncharacteristic, so u must vanish in a neighborhood 
of Xr i£ ; hence R — D. The regions X r<6 exhaust X, completing the proof. 

□ 

A proof of this lemma can also be found in [KaKuLa] . 

Wc finish the proof of Theorem 12. li as follows: The function c(x) is bounded away 
from 0, and all of its first and second derivatives are bounded functions. This implies 
gij is a metric whose curvature components are bounded, and as a result, there exists 
a number 5 > such that Vx G K™ there are geodesic normal coordinates about x 
covering a region {y | d(x, y) < 5}. We obtain the following corollary of Lemma T3. 21 

Corollary 3.3. Let r > 0. If u = on B(p,r) x [—H,H], then u — on the set 



(32) Y = {(x,i) | d(x,z) <r + vam(H-\t\,S)} . 

Proof. Let (x,£) £ Y. In the case when t G (—H,—H + S), wc apply Lemma [3.21 
knowing u vanishes on B(p,r) x \—H, —H + 25] to conclude u(x,t) = 0. Similarly, 
if t G (H — 5, H) we apply the lemma knowing u vanishes on B(p, r) x [H — 26, H]. 
Finally, if —H + 5 < t < H — 8, wc find u(x, t) = by applying the lemma, knowing 
u vanishes on B(p, r) x [t — 5, t + 6} . □ 

Repeated applications of this corollary to the cylinders 



B{p, p + nS)x [-(H - nS), H - n5] 

show that if u vanishes on B(p, p) x [-H, H] for some p > 0, then u vanishes at (x, t) 
if d(x, z) + \t\ < H . In particular, f(xo) = u(xq, 0) = 0. Hence / = 0, concluding the 
proof of Theorem l2.ll □ 

Remark 3.4. This work was done several months ago and was included my 
preliminary Ph.D. examination in June, 2008. It appears my results are similar to 
the results contained in a paper of P. Stefanov and G. Uhlmann 
( Utp://arxiv.org/PS-cache/arxiv/pdf/0902/0902.1973vl.pdf$ 
that was just announced. 
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